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Abstract 

Background: Emerging whitefly transmitted begomoviruses are major pathogens of vegetable and fibre crops 
throughout the world, particularly in tropical and sub-tropical regions. Mutation, pseudorecombination and 
recombination are driving forces for the emergence and evolution of new crop-infecting begomoviruses. Leaf curl 
disease of field grown radish plants was noticed in Varanasi and Pataudi region of northern India. We have 
identified and characterized two distinct monopartite begomoviruses and associated beta satellite DNA causing 
leaf curl disease of radish {Rophonus sotivus) in India. 

Results: We demonstrate that RaLCD is caused by a complex of two Old World begomoviruses and their 
associated betasatellites. Radish leaf curl virus-Varanasi is identified as a new recombinant species, Radish leaf 
curl virus (RaLCV) sharing maximum nucleotide identity of 87.7% with Tomato leaf curl Bangladesh virus- 
[Bangladesh:2] (Accession number AF188481) while the virus causing radish leaf curl disease-Pataudi is an 
isolate of Croton yellow vein mosaic virus-[lndia] (CYVMV-IN) (Accession number AJ507777) sharing 95.8% 
nucleotide identity. Further, RDP analysis revealed that the RaLCV has a hybrid genome, a putative 
recombinant between Euphorbia leaf curl virus and Papaya leaf curl virus. Cloned DNA of either RaLCV or 
CYVMV induced mild leaf curl symptoms in radish plants. However, when these clones (RaLCV or CYVMV) were 
individually co-inoculated with their associated cloned DNA betasatellite, symptom severity and viral DNA 
levels were increased in radish plants and induced typical RaLCD symptoms. To further extend these studies, 
we carried out an investigation of the interaction of these radish-infecting begomoviruses and their associated 
satellite, with two tomato infecting begomoviruses {Tomato leaf curl Gujarat virus and Tomato leaf curl New 
Delhi virus). Both of the tomato-infecting begomoviruses showed a contrasting and differential interaction with 
DNA satellites, not only in the capacity to interact with these molecules but also in the modulation of 
symptom phenotypes by the satellites. 

Conclusion: This is the first report and experimental demonstration of Koch's postulate for begomoviruses 
associated with radish leaf curl disease. Further observations also provide direct evidence of lateral movement of 
weed infecting begomovirus in the cultivated crops and the present study also suggests that the exchange of 
betasatellites with other begomoviruses would create a new disease complex posing a serious threat to crop 
production. 
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Background 

Geminiviruses (family Geminiviridae) have circular sin- 
gle stranded DNA genomes that are encapsidated in 
twinned quasi-isometric particles and are classified into 
four genera: Begomovirus, Mastrevirus, Curtovirus and 
Topocuvirus [1,2]. The largest genus, Begomovirus com- 
prise of viruses transmitted by whitefly (Bemisia tabaci 
Genn.). Begomoviruses are assumed to have been co- 
evolving with their hosts for a long time, however, it is 
in the past two decades, these viruses have become an 
economically important plant pathogens [3-5]. These 
viruses contain genomes consisting of either one or two 
similar-size DNA components [6]. The DNA- A compo- 
nent encodes all viral functions required for replication, 
control of gene expression, encapsidation and vector 
transmission [7,8]. The DNA-B, code for two proteins 
which are involved in movement of the virus between 
and within plant cells [6] . Both components of bipartite 
begomoviruses are required for systemic infection and 
symptom induction [9]. In contrast, monopartite bego- 
moviruses such as isolates of Tomato yellow leaf curl 
virus, Tomato leaf curl virus and Cotton leaf curl Mul- 
tan virus [10-13] possess only a single genomic compo- 
nent resembling DNA-A, which alone is capable of 
inducing disease symptoms. However, for some other 
monopartite begomoviruses like Ageratum yellow vein 
virus, Bhendi yellow vein mosaic virus, Cotton leaf curl 
Multan virus, Eupatorium yellow vein virus, Tomato yel- 
low leaf curl China virus and Cotton leaf curl Gezira 
virus, association of betasatellite (DNA-P) has been 
found to be essential for the induction of typical disease 
symptoms [14-19]. Betasatellite is a circular, single- 
stranded DNA molecule of -1.35 kb length with a single 
open-reading frame (ORF) PCI, an adenine-rich region, 
a satellite conserved region having nonanulceotides 
(TAATATTAC) and it shares negligible sequence simi- 
larity with either DNA-A or DNA-B of bipartite bego- 
moviruses [20-22]. 

Radish (Raphanus sativus L.; Yamily-Brassicaceae) is 
an important vegetable crop grown throughout India. 
Leaf curl disease of radish (RaLCD), for the first time, 
was reported from Punjab province of Pakistan [23]. In 
India, RaLCD was first reported from a homestead gar- 
den and adjoining farmer's field near Varanasi [24]. 
Later, RaLCD was also observed in Pataudi, Haryana, 
India. The infected plants remain stunted and leaves 
exhibited both upward and downward curling along 
with conspicuous vein enations. 

In the present study, we have characterized two dis- 
tinct begomoviruses and betasatellites associated with 
plants showing symptoms of RaLCD in Varanasi and 
Pataudi, India. These begomovirus species, together with 
the betasatellite can induce RaLCD. We further 



demonstrate the complexity of the interactions between 
RaLCD-associated satellites with radish and tomato- 
infecting begomoviruses indicating co-adaptation of 
weed infecting begomovirus in cultivated crop species 
like radish, hitherto a non-host. 

Results 

Cloning and genome organization of two begomovirus 
species associated with RaLCD in India 

RaLCD symptoms developed on grafted Raphanus sati- 
vus plants were identical to those of naturally infected 
plants. The initial symptoms appeared as downward cur- 
ling on the young leaves about 3-4 weeks after grafting, 
later infected leaves exhibited typical upward and down- 
ward leaf curling, enation on adaxial side and stunted 
growth of plants, which also failed to bear any flower. 
Southern blot hybridization and PCR based detection 
using begomoviruses specific primers [25] revealed asso- 
ciation of begomovirus with RaLCD infected plants from 
Varanasi, Uttar Pradesh and Pataudi, Haryana in India. 
Full-length (-2.7 kb) viral genome from both samples 
(Cholapur and Pataudi) were cloned. Several putative 
full-length clones were obtained and preliminary restric- 
tion analysis indicated the presence of single class of 
molecule in both samples (data not shown). Two full- 
length molecules from each sample were selected and 
sequenced. The complete sequence of RaLCD-asso- 
ciated- Varanasi virus (GenBank EF 175733) and RaLCD- 
associated-Pataudi virus (GenBank FJ593629) isolate 
were of 2756 nucleotides (nt) and 2759 nt, respectively. 
The genome organization of both the isolates was similar 
to that of Old World monopartite begomoviruses [i.e. 
two virus-sense ORFs (VI and V2) and four complemen- 
tary-sense ORFs (CI, C2, C3, C4)]. Intergenic region (IR) 
sequences (-290 nt) contained conserved nonanucleotide 
sequence, putative Rep (CI) protein binding sites and the 
TATA box. The nucleotide identity of full-length genome 
and IR sequences between RaLCD-associated-Varanasi 
and RaLCD-associated-Pataudi were 79.6 and 71.6%, 
respectively, indicating involvement of two distinct spe- 
cies (Table 1). 

The sequence of the begomovirus DNA component of 
the RaLCD-associated-Varanasi virus had maximum 
nucleotide identity with Tomato leaf curl Bangladesh 
virus- [Bangladeshi] (ToLCBDV-[BD:2]) (-87.7%) (Table 
1). Phylogenetic analysis also support that RaLCD-asso- 
ciated-Varanasi is a distinct species, having very low 
bootstrap values with the Euphorbia leaf curl virus- 
[China:Guangxi 35:2002] (EuLCV-[CN:Gx35:02]) (Figure 
1). Thus, RaLCD-associated-Varanasi virus is an isolate 
of a new begomovirus species, and the name Radish leaf 
curl virus- [India: Varanasi:2005] (RaLCV-[IN:Var:05]) is 
given following (ICTV) guidelines [2]. 
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Table 1 Nucleotide identities (%) for full-length genome 
(nt) and IR and amino acid (aa) sequence identities (%) 
for ORFs of two begomoviruses associated with RaLCD in 
India and selected previously characterized 
begomoviruses 



Begomoviruses 


DNA- 
A 


Rep 


TrAP 


REn 


AC4 


CP 


AV2 


ICR 


RaLCV-[IN:Var:05] 


1 0LLdDV-|_dU:zJ 


87.7 


84.8 


80.6 


86.6 


54.6 


94.1 


98.8 


75.3 


TbCSV-[CN:Yn35:01] 


87.0 


84.5 


87.3 


88.1 


50.5 


93.9 


99.2 


71.2 


AEV-[NP:01] 


86.9 


84.5 


85.8 


88.1 


65.9 


93.9 


98.4 


74.0 


TbCSV-[CN:Yn282: 
Age:03] 


86.9 


83.9 


88.1 


88.1 


48.0 


93.0 


99.2 


71.2 


PaLCuV-IN[IN:Luc] 


86.4 


73.6 


96.3 


97.0 


44.7 


91.5 


97.7 


58.0 


CLCuKV-Man[PK: 
M806b:96] 


85.8 


87.5 


84.3 


82.8 


66.0 


96.6 


98.8 


40.2 


EuLCV-[CN:Gx35:02] 


85.7 


90.0 


82.1 


88.1 


86.5 


78.0 


93.8 


69.7 


ChiLCV-IN[IN::05] 


84.0 


86.4 


82.1 


73.1 


52.9 


98.8 


89.3 


75.3 


ChiLCV-Mul[PK: 
Mul:98] 


83.0 


84.5 


79.4 


80.6 


53.6 


87.3 


96.1 


64.4 


ToLCGV-[IN:Var:01] 


80.7 


83.1 


79.9 


77.6 


51.0 


86.1 


80.9 


68.4 


CYVMV-[IN] 


80.4 


75.9 


83.6 


87.3 


48.2 


87.3 


83.2 


59.1 


BYVMV-IN[IN:Mad] 


74.9 


81.2 


64.9 


67.9 


69.6 


70.3 


93.4 


38.7 


ToLCNDV-IN[IN:ND: 
Svr:92] 


74.2 


77.8 


60.4 


66.4 


62.1 


73.2 


93.0 


60.0 


SLCCNV-IN[IN:Luc: 
Pum] 


73.6 


75.6 


57.5 


64.8 


53.5 


71.4 


93.0 


40.3 


MYMIV-[IN:ND: 
Bg3:91] 


63.5 


70.6 


52.2 


40.2 


45.5 


50.0 


72.2 


10.9 


CYVMV-[IN:Pat: 
Rad:08] 


RaLCV 


79.6 


75.6 


82.8 


89.6 


44.7 


83.2 


87.3 


71.6 


ToLCBDV-[BD:2] 


80.1 


77.0 


85.8 


84.3 


54.6 


82.8 


91.5 


67.5 


TbCSV-[CN:Yn35:01] 


80.0 


76.7 


77.6 


85.1 


35.3 


83.2 


92.2 


71.3 


AEV-[NP:01] 


80.9 


77.8 


79.9 


86.6 


52.9 


82.4 


91.3 


70.5 


TbCSV-[CN:Yn282: 
Age:03] 


79.9 


76.2 


78.4 


85.1 


34.1 


83.2 


91.3 


71.3 


PaLCuV-IN[IN:Luc] 


86.0 


92.2 


81.3 


88.1 


95.3 


82.4 


97.7 


59.7 


CLCuKV-Man[PK: 
M806b:96] 


78.1 


75.0 


86.6 


87.3 


44.7 


83.6 


98.8 


54.5 


EuLCV-[CN:Gx35:02] 


78.1 


77.7 


85.8 


88.8 


45.9 


79.7 


72.9 


73.3 


ChiLCV-IN[IN::05] 


78.4 


77.3 


82.1 


77.6 


35.3 


82.4 


89.0 


70.1 


ChiLCV-Mul[PK: 
Mul:98] 


78.3 


76.5 


82.4 


84.3 


31.8 


82.0 


83.1 


74.6 


ToLCGV-[IN:Var:01] 


80.5 


76.5 


83.6 


81.3 


35.3 


96.5 


85.2 


66.9 


GYVMV-PN] 


95.8 


94.7 


96.3 


93.3 


90.6 


100 


98.3 


94.4 


BYVMV-IN[IN:Mad] 


71.3 


70.0 


64.2 


71.6 


40.0 


80.5 


66.9 


54.2 


ToLCNDV-IN[IN:ND: 
Svr:92] 


71.5 


69.5 


64.9 


64.2 


34.5 


82.8 


68.8 


64.2 


SLCCNV-IN[IN:Var: 
Pum] 


70.8 


67.3 


60.4 


68.0 


37.9 


81.2 


67.0 


61.5 


MYMIV-[IN:ND: 
Bg3:91] 


62.0 


65.8 


50.0 


40.2 


30.6 


74.1 


47.3 


46.2 



"Identity over 80% is Bold 



However, RaLCD-associated-Pataudi virus sequence 
was 95.8% identical to Croton yellow vein mosaic virus- 
[India] (CYVMV-[IN]) isolates (Table 1) and therefore, 
it can be considered as a variant of CYVMV-[IN] and 
designated as Croton yellow vein mosaic virus- [India: 
Pataudi:Radish:2008]. 

The predicted amino acid sequences of the ORFs of 
RaLCV- [IN: Var:05] were compared to those of other 
begomoviruses. TrAP, REn, CP and AV2 of RaLCV- [IN: 
Var:05] showed maximum amino acid identity (91-98%; 
Table 1) with Papaya leaf curl virus-India [India:Luck- 
now] (PaLCuV-IN[IN:Luc]) with which full-length 
DNA-A shared 86.4% identity (Table 1). However, AC4 
and Rep shared maximum amino acid (aa) identity 
(86.5% and 90%, respectively) (Table 1) with weed- 
infecting begomovirus, EuLCV-[CN:Gx35:02]. Interest- 
ingly, IR of DNA-A showed maximum identity with 
ToLCBDV-[BD:2] (75.3%) (Table 1). Such divergence in 
the amino acid sequence of the encoded proteins pre- 
sent strong evidence of recombination in the genome of 
RaLCV- [IN: Var:05]. Sequences of the CYVMV-[IN:Pat: 
Rad:08] VI (CP), V2, CI (Rep), C2, C3 and IR ORFs 
were generally > 95% [nt and aa] identical with 
CYVMV-[IN], whereas maximum identity of C4 
sequences (90% nt and 85% aa) were with PaLCuV-IN 
[IN:Luc] (Table 1). 

A multiple alignment of the IR of six begomoviruses 
revealed that RaLCV- [IN: Var:05] and ToLCBDV-[BD:2] 
share very related iterons sequence (GGTGA-AC/T- 
GGTAC) whereas CYVMV-[IN:Rad:Pat:08] and 
CYVMV-IN iterons sequences were identical (GGGGA- 
CTC-GGGGGA) (data not shown). The major difference 
between the RaLCV- [IN: Var:05] and ToLCBDV-[BD:2] 
was found in spacer sequence of two iterons, as RaLCV- 
[IN:Var:05] contains AC while ToLCBDV-[BD:2] con- 
tains AT. A third iteron was present in 5' end of the 
intergenic common region (ICR) of RaLCV. Noticeably, 
besides GC rich region, stem loop region and the TATA 
Box were conserved for all geminiviruses and all the 
viruses tested here shared two short sequences, 
AATGGCA and TAAT, located at the 5'end of the 
region between TATA box and stem-loop. 

Cloning of betasatellite DNAs and genome organization 

Despite repeated attempts, we could not identify any 
DNA-B molecule from the infected samples, indicating 
monopartite nature of virus isolates. Full-length genome 
(-1.3 kb) of betasatellites was cloned from both (Vara- 
nasi and Pataudi) RaLCD samples using universal pri- 
mers [26]. The complete nucleotide sequence of 
betasatellites associated with RaLCV- [IN: Var:05] and 
CYVMV-[IN:Pat:Rad:08] was determined as 1358 nt 
(GenBank Accession no EF175734) and 1367 nt (Gen- 
Bank Accession no. FJ593630), respectively. These two 
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PaLCuV-IN[IN:Luc] ;Y15934 

PaLCuV-PK[PK:Cot:02] ;AJ436992 

RaLCV-[IN:Var:03] ;EF175733 

EuLCV-[CN:Gx35:02] ;AJ553121 

ToLCBDV-[BD:2] ;AF1 38481 

AEV-[NP:01] ;AJ437618 

r TbCSV-[CN:Yn35:01] ;AJ420318 
Wl; TbCSV-[CN:Yn282:Age:03] ;AJ971266 

ToLCGV-[IN:Var:01] ;AY1 90290 

Tol_CKV-Ban[IN:Ban:93] ;U38239 

CLCuKV-Man[PK:M806b:96] :AJ002449 

CLCuKV-Man[IN:Dab];AY456683 

BYVMV-IN[IN:Mad] ;AF241479 

OYVMV-[PK:Fai201:95] ;AJ002451 

CLCuMV-Fai[PK:Y62:95] ;AJ002447 

CLCuMV-IN[IN:Sri:94];AF36301 1 

Tol_CSLV-[LK:Ban:97] ;AF274349 

ToLCBV-A[IN:Ban1] ;Z48182 

Tol_CBV-[IN:Pun:05] AY7 54814 

SLCCNV-IN[IN:Luc:Pum] ;DG026296 

ToLCNDV-IN[IN: ND:Svr:92] ; U1 501 5 

MYMIV-[IN:ND:Bg3:91] ;AF126406 



—I 

01 



0.15 



0.10 



0.05 



Nucleotide Substitution (xlOO) 
Bootstrap Trials = 1000, seed = 1 1 1 

Figure 1 Phylogenetic relationship of full-length DNA-A of radish-infecting begomoviruses with other begomoviruses. Phylogenetic 
tree was generetaed based on aligned complete DNA-A sequences using ClustalW. All the sequences were obtained from GenBank. The viruses 
understudy are indicated as bold. 



betasatellites had structural features similar to those of nucleotide sequence identity between RaLCV-[IN: 
other betasatellites i.e. a single ORF pel, satellite con- Var:05] and CYVMV-[IN:Pat:Rad:08] was 63.2%, indicat- 
served region (SCR) and A rich region. The complete ing that these are distinct betasatellite species. The 
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complete sequence of DNA betasatellite associated with 
RaLCV-[IN:Var:05] and CYVMV-[IN:Pat:Rad:08] shared 
maximum nucleotide identity with Tobacco leaf curl 
betasatellite [PakistamRahim Yar Khan: 1998] (TbLCB- 
[PK:RYK:98]) (83%) and Croton yellow vein mosaic 
betasatellite [Pakistan:Punjab:2006] (CroYVMB-[PK: 
Pun:06]) (91.1%), respectively (Table 2). Phylogenetically, 
RaLCV-[IN:Var:05] associated DNA betasatellite clusterd 
with the one associated with tobacco disease with 
TbLCB-[PK:RYK:98], whereas the CYVMV-[IN:Pat: 
Rad:08] associated DNA betasatellite clustered with 
another isolate of this species. (Figure 2). 

ORF (3C1 of RaLCB-[IN:Var:05] and CYVMB-[IN:Pat: 
Rad:08] showed maximum nucleotide identity with 
TbLCB-Pak-[NIB 12-1] (89.1%) and CroYVMB-[PK: 
Pun:06] (93.8%), respectively (Table 2). Similarly, ORF 
PCI encoded protein of RaLCV-[IN:Var:05] and 



Table 2 Per cent identity (nucleotide) among DNA-(3 of 
two Radish leaf curl virus isolates and with selected 
begomoviruses originating in Asia 



Beta 


Percent 


CI 


SCR 


TbLCB- Radish 


TbLCB-[PK:RYK:98] 


83.0 


89.1(81.1) 


84.4 


ToLCB-[IN:Var:06] 


72.7 


79.3(76.3) 


93.2 


TYLCTHB-[IN:Aur:06] 


72.3 


76.6(75.0) 


93.0 


AYVSLB-[IN:Mad:03] 


68.2 


75.9(71.2) 


90.0 


TbCSB-[CN:Yn 115:02] 


67.2 


68.2(68.6) 


93.9 


ToYLCCNB-[CN:Yn261:02] 


66.2 


69.6(61.9) 


90.8 


AYLCuB-[PK:Fai3:94] 


66.0 


72.6(56.5) 


94.3 


CroYVMB-[PK:Pun:06] 


62.5 


62.6(57.5) 


91.7 


ToLCB-[IN:ND:02 


62.0 


72.1(66.9) 


80.3 


PaLCuB-[IN:Chi:05] 


61.3 


64.5(59.8) 


85.5 


PaLCuB-[IN:Jab:03] 


61.2 


65.7(61.9) 


64.4 


ToLCB-[IN:CP:04] 


61.1 


70.8(66.9) 


84.1 


CYVMB-Radish 


RaLCB-Var 


63.2 


64.7 (61.9) 


90.8 


TbLCB-[PK:RYK:98] 


60.4 


65.0 (61.0) 


83.8 


ToLCB-[IN:Var:06] 


62.4 


64.7 (57.6) 


91.9 


TYLCTHB-[IN:Aur:06] 


61.8 


65.0 (56.8) 


90.8 


AYVSLB-[IN:Mad:03] 


59.7 


66.7 (60.2) 


88.2 


TbCSB-[CN:Yn 115:02] 


59.8 


63.3 (57.6) 


91.2 


ToYLCCNB-[CN:Yn261:02] 


60.0 


64.7 (60.2) 


86.4 


AYLCuB-[PK:Fai3:94] 


62.3 


66.0 (62.3) 


90.8 


CroYVMB-[PK:Pun:06] 


91.1 


93.8 (90.7) 


95.6 


ToLCB-[IN:ND:02] 


59.8 


65.8 (61.0) 


85.9 


PaLCuB-[IN:Chi:05] 


60.8 


70.9 (66.9) 


87.6 


PaLCuB-[IN:Jab:03] 


61.3 


71.4 (66.1) 


84.1 


ToLCB-[IN:CP:04] 


59.8 


63.9 (61.0) 


85.8 



Figures in bracket indicate putative deduced amino acid sequence identity. 
"Identity over 80% is Bold 



CroYVMB-[IN:Pat:Rad:08] associated DNA betasatellite 
showed maximum identity with TbLCB-Pak-[NIB 12- 
1] and CroYVMB-[PK:Pun:06] (87.3% and 90.7%, 
respectively) (Table 2). The conserved ORF PCI of 
RaLCB-[IN:Var:05] and CroYVMB-[IN:Pat:Rad:08] 
encoded a 13.7 kDa protein. The nucleotide sequences 
of SCR of RaLCV-[IN:Var:05] and CYVMV-[IN:Pat: 
Rad:08] associated DNA betasatellite showed > 90% 
identity with those of other DNA betasatellites (Table 
2). SCR of RaLCV-[IN:Var:05] associated DNA betasa- 
tellite and CYVMV-[IN:Pat:Rad:08] associated DNA 
betasatellite showed maximum identity with AYLCuB 
(94.3%) and CroYVMB-[PK:Pun:06] (95.6%), respec- 
tively (Table 2). 

Detection of recombination 

RaLCV-[IN:Var:05] (RLA) and CYVMV-[IN:Pat:Rad:08] 
(CRA) genomes were analyzed using two independent 
methods (simplot graph, and RDP analysis) and the 
results indicated definite evidence of recombination 
within RaLCV-[IN:Var:05] genome (Figure 3). In 
RaLCV-[IN:Var:05], two significant breakpoints detected 
at nts 725 and 1549 with significant 'P value' (0.01) (Fig- 
ure 3c), corresponding to the region of the genome with 
> 90% identity with the sequence of PaLCuV. The N- 
terminal amino acid sequence of Rep was most similar 
to EuLCV (87%) than other geminiviruses. Evidence of 
recombination with EuLCV like virus was found within 
Rep coding region (from nt 2154 to 2519). Bootscan and 
Simplot analysis also supported these observations (Fig- 
ure 3a, b). To confirm recombinant nature of RaLCV- 
[IN:Var:05] genome, dendrograms were also generated 
based on nucleotide sequence of virion sense ORF 
(AVI), complimentary sense ORFs (AC1, AC2 and AC3) 
with that of related begomoviruses. Interestingly, 
RaLCV-[IN:Var:05] clustered with PaLCuV when AVI, 
AC2 and AC3 sequences were used for phylogenetic 
analysis, whereas RaLCV-[IN:Var:05] was placed with 
EuLCV in the same clade in the dendrogram generated 
using AC1 nucleotide sequences (data not shown). 

To detect possible evidence of recombination in the 
CYVMV-[IN:Pat:Rad:08] genome, Simplot and RDP ana- 
lysis were performed with begomoviruse bearing the clo- 
sest identity viz., CYVMV-[IN], and with two tomato- 
infecting begomoviruses viz., Tomato leaf curl Gujarat 
virus- [India: Varanasi:2001] (ToLCGV-[IN:Var:01]) and 
Tomato leaf curl New Delhi virus-India [India:New 
Delhi:Severe:1992] (ToLCNDV- [IN:ND:Svr:92]). Sim- 
plot and RDP analyses revealed that this is not a recom- 
binant virus; rather it is a derivative of CYVMV-[IN]. 
Similarly, simplot and RDP analysis carried out for beta- 
satellite molecules DNAs (RaLCB-[IN:Var:05] and 
CroYVMB-[IN:Pat:Rad:08]) could not reveal any evi- 
dence of recombination. 
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Figure 2 Phylogenetic relationship of RaLCB and CroYVMB with betasatellites associated with other begomoviruses. Phylogenetic tree 
was geenerated based on alignment of complete genome sequences using ClustalW. The viruses understudy are indicated as bold. 



Infectivity of cloned DNAs 

Partial tandem repeats of DNA-A and DNA betasatellite 
were constructed in binary vector and their infectivity 
were tested (either alone or together in combination) 



onto radish, tomato and N. benthamiana plants by 
agroinoculation. Radish plants co-inoculated with RLA 
+ RLp developed systemic symptoms including down- 
ward leaf curling within 20 days post inoculation (dpi) 
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(Figure 4a, b; Table 3), however, RLA inoculated plants 
exhibited milder and delayed (6 days) symptoms as 
compared to infection by RLA + RLp (Figure 4a; Table 

3) . Similarly, radish plants co-inocluated with CRA + 
CRp exhibited severe leaf curling, twisting of petiole and 
enations symptoms (Figure 4a, d; Table 3), while CRA 
inoculated plants showed delayed (7 days) and milder 
symptom (Figure 4Ac; Table 3). The symptoms observed 
under experimental studies on radish plants were very 
similar to those observed under field conditions. The 
presence of viral genomic DNA and DNA betasatellites 
in infected radish plants was confirmed by Southern 
blot (Figure 4b and 4c) and PCR analysis (results not 
shown). 

Tomato plants agroinoculated with the RLA and RLp 
did not develop any symptom (Figure 4a, b), whereas 
those co-inoculated with CRA and CRp developed 
severe leaf curling and stunting symptoms. CRA infected 
plants produced systemic leaf curling on upper leaves 
and mild stunting at later stages (Figure 4;able 4), but 
when plants were co-inoculated with CRA and CRp 
they were severely stunted (Figure 4; Table 4). Agroino- 
culation of infectious clones of RLA and CRA to N. 
benthamiana resulted in leaf curling, enations and mild 
stunting of test plants (Figure 4Aal and 4Acl) followed 
by typical severe leaf curling, twisting of petiole, interve- 
inal chlorosis, severe stunting and enations on underside 
of leaves when co-inoculated with RLp and CRp, respec- 
tively (Figure 4Abl and 4Adl; Table 4). Presence of viral 
genomic DNA and DNAps in infected N. benthamiana 
and tomato were confirmed by Southern blot (Figures 
4D, E, 5a and 5b) and PCR analysis (results not shown). 
Accumulation of DNA component (RLA and CRA) was 
enhanced when DNA betastaellite was associated. 

Studies to determine whether RLp is trans-replicated 
by CRA and vice versa, agroinoculation experiments 
were conducted. Results suggested that RLp could be 
transreplicated by CRA and plants co-infected with CRA 
and RLp develop severe leaf curling, leaf crinkling, inter- 
veinal chlorosis, enation on underside of leaves and 
severe stunting (Table 4). Similarly, RLA could also 
trans-replicate with CRp and this combination produce 
distinct symptoms in N. benthamiana. Therefore, both 
the DNAps could be replicated by each other helper 
begomovirus leading to severe symptom phenotype 
(Table 4). 

RLA and CRA do not transreplicate DNA-B component of 
ToLCNDV and ToLCGV 

Co-inoculation of RLA and CRA with either DNA-B of 
ToLCNDV (NB) or ToLCGV (VB) did not influence 
symptom severity nor the time required for symptom 
appearance (Figure 6a-al, 6b-bl, f-fl and 6g-gl; Table 

4) . Thus, the symptom expression pattern was similar as 



observed in plants inoculated with RLA and CRA alone. 
In addition, presence of DNA-B could not be ascer- 
tained for both the cases (Figure 4D-E and 5d). It is 
relevant to mention here that iteron sequences are dif- 
ferent among RLA, CRA with both the DNA-B compo- 
nents of tomato-infecting begomoviruses. This could be 
a barrier for efficient transreplication of the heterolo- 
gous DNA-B component. 

RLp and CRp can be a substitute for DNA B of either 
ToLCNDV or ToLCGV 

ToLCNDV is a bipartite begomovirus that requires both 
DNA components for systemic infection [27]. Tomato 
plants inoculated with DNA-A (NA) and DNA-B (NB) 
of ToLCNDV developed leaf curling within 8 dpi and 
later the infected plants showed stunted growth along 
with severe leaf curling (30 dpi) (Figure 6ql), whereas 
tomato plants inoculated with only DNA-A remained 
symptomless (Table 4). Plants inoculated with NA alone 
lacked detectable levels of DNA in the newly emerged 
leaves when analyzed by PCR and Southern blot at 30 
dpi (data not shown). To test whether a DNA betasatel- 
lite could influence the pattern of NA infection, tomato 
plants were inoculated with a mixture of NA + RLp or 
NA + CRp. A proportion (58%) of plants inoculated 
with NA + RLp showed initial leaf curling (18 dpi) that 
later turned into pronounced leaf curling and enations 
at 30 dpi (Figure 6dl, Table 4). Southern blot hybridiza- 
tion confirmed the presence of replicative forms of both 
NA and RLp in the distal leaves of symptomatic plants 
(Figure 4E lane 6 and G). However, accumulation of NA 
was lower in these plants as compared to plants inocu- 
lated with NA + NB (Figure 4G). Interestingly, NA + 
CRP inoculated tomato plants remain symptomless and 
newly emerging leaves lacked detectable levels of either 
NA or CRp (Figures 5b-c and 6il.). N. benthamiana 
plants inoculated with NA + RLp showed distinguish- 
able symptoms like typical downward leaf curling, inter- 
veinal chlorosis and twisting of petiole (Figure 6d), 
whereas NA + CRp induced leaf curling, twisting of 
petiole and leaf puckering (Figure 6i). Accumulation of 
viral DNA (NA) was less (-51%) in plants inoculated 
with NA + RLp or NA + CRp as compared to those 
inoculated with homologous combination (NA + NB) 
(Figures 4F and 5c). Hence, RLp and CRp could support 
efficient movement of NA in N. benthamiana. However, 
CRP cannot facilitate efficient movement of NA on 
tomato unlike RLp. 

ToLCGV is a mono-bipartite species [28,29] and 
DNA-A (VA) alone is capable of causing systemic infec- 
tion on N. benthamiana and tomato, although associa- 
tion of DNA-B results in increased symptom severity 
[28,29]. Neither RLp nor CRp could enhance systemic 
accumulation of VA as compared to VA + VB (Figures 
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Figure 4 Symptom expression and DNA accumulation in test plants. Plant phenotypes [A] radish, N. benthomiano and tomato seedling 
following inoculation with a RaLCV DNA- A (RLA) and CYVMV-[IN:Pat:Rad:08] DNA- A (CRA)were inoculated either alone or with RaLCV DNA- 
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relative accumulation of viral DNA in inoculated [B-C] radish, [D] N. benthamiana (left panel)and [E] tomato (right panel) plants. 



Singh et al. Virology Journal 201 2, 9:43 Page 1 0 of 1 8 

http://www.virologyj.eom/content/9/1/43 



Table 3 Infectivity of viral genome and DNA betasatellites on radish plants 

Test plant Construct combination used for Plants showed symptom/Plants Symptom First symptom appearance 

(s) agroinoculation inoculated (dpi) 

Radish RLA 11/20 MLC 26 

RLA + RLJ3 15/20 SL, LC, Ens 20 

CRA 12/20 MLQLR 24 

CRA+CRJ3 1 7/20 LC, St, Ens, 1 7 

Tw 

VA+VB 0/16 

NA+NB 0/16 

RLA, RaLCV; RLp, RaLCB; CRA, CYVMV; CRp, CroYVMB; NA, ToLCNDV DNA-A; NB, ToLCNDV DNA-B; VA ToLCGV-DNA-A; VB, ToLCGV-DNA-B; MLC mild leaf curling; SL 
small leaves; LC leaf curling; Ens enation; LR Leaf rolling; Tw twisting of petiole 



4F-G, and 5e), although inoculated plants exhibited 
severe symptom (Table 4). For example, in N. benthami- 
ana, presence of either RLp or CRp alongwith VA 
together caused twisting of the petiole, interveinal 
chlorosis and more stunting of the plants, while tomato 
plants exhibited stunting alongwith severe leaf curling as 
compared to plants inoculated with VA + VB or VA + 
CRP (Figure 6c-cl, 6h-hl and 6p-pl, Table 4). 

Synergism between CRA, ToLCNDV and ToLCGV 

Since, CRA + CRp could infect tomato and produced 
distinct symptoms like leaf curling, enation etc. (Table 
4), we tested the potential consequences/interaction of 
these molecules with two predominant tomato-infecting 
begomoviruses (ToLCNDV and ToLCGV) from north- 
ern India. For this, tomato and N. benthamiana plants 
were inoculated with various combinations of viral DNA 
(Table 4). Plants inoculated with heterologous combina- 
tions of begomovirus species showed unusually severe 
symptoms as compared with plants infected with homo- 
logous combinations(CRA + CRp or VA + VB or NA + 
NB) (Figure 6e-q). More interestingly, symptom severity 
was enhanced whenever plants were inoculated with 
CRA + CRp alongwith either NA or VA. For example, 
when CRA or CRA + CRp were mixed with VA or NA 
or VA + VB or NA + NB, plants exhibited severe symp- 
toms (Figure 6j-o). Plants inoculated with both compo- 
nents of any two virus species exhibited symptoms at 
least a day earlier than those inoculated with CRA + 
CRp (Table 4). To investigate synergism at the level of 
viral DNA accumulation, Southern blot analysis was 
conducted with probes specific for each of the six DNA 
components (Figure 5). The level of CRA in plants 
inoculated with CRA + VA + CRp, CRA + NA + CRp, 
CRA + VA + VB and CRA + NA + NB was low in 
tomato as compared to N benthamiana (Figure 5a). 
However, the level of CRA was -1.4 times higher when 
VA + VB was also inoculated along with CRA + CRp 
(Figure 5a), indicating CRp assisted synergistic interac- 
tion. Unlike VA where no increase was observed, NA 



accumulation was marginally increased in plant inocu- 
lated with all four molecules (Figure 5c and 5e). Inter- 
estingly, level of either VB or NB remained unchanged 
in plants inoculated with any of the following combina- 
tions CRA + VA + VB, CRA + NA + NB, CRA + VA + 
VB + CRp, CRA + NA + NB + CRp as compared to 
either VA + VB or NA + NB, respectively (Figure 5d 
and 5f). The level of CRp was reduced by 53% and 42% 
when it was co-inoculated with VA and NA, respectively 
as compared with the wild-type combination (Figure 
5b). Accumulation of CRp decreased upto 50-60% in 
plants inoculated with either CRA + VA + VB + CRp or 
CRA + NA + NB + CRp, in comparison to plants inocu- 
lated with CRA + CRP alone (Figure 5b). Together, 
these results indicate that CRA may interact synergisti- 
cally with VA + VB in the presence of CRp. 

Discussion 

In India, radish was not considered as a host for bego- 
moviruses till Singh et al. [24] reported association of a 
begomovirus with leaf curl disease from Varanasi region. 
Subsequently, RaLCD was also observed from Pataudi, a 
place located -800 km West from Varanasi. This indi- 
cates emergence of RaLCD as a new disease in India, 
incidence of which was recorded between 10-40%. Here, 
we ascertain that RaLCD is caused by a complex of two 
Old World monopartite begomovirus species and two 
distinct betasatellite molecules. On the basis of genome 
organization and phylogenetic analysis, the virus from 
Varanasi region (RaLCV- [IN: Var:05]) is considered as a 
novel begomovirus species, which share maximum 
sequence identity 87.7% with ToLCBV-[BD], whereas 
virus from Pataudi share maximum sequence identity 
95% with CYVMV-IN and is considered as a CYVMV- 
IN variant on the basis of DNA-A sequence identity cri- 
teria of ICTV [2]. Generally, begomovirus isolates shar- 
ing less than 89% nucleotide sequence identity are 
considered to be distinct species, sharing 90-92% are 
considered as new strain and those with more than 94% 
identity and are considered to be variant of that strain 
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Table 4 Constructs, their infectivity and symptoms induced on agroinoculated plants 



Test plant 

(s) 


Construct combination used for 
agroinoculation 


Plants showed symptom/Plants 
inoculated 


Symptom 


First symptom 
appearance (dpi) 


N. 

benthamiana 


RLA 


7/16 


LC, Ens 


20 




RLA + Rl_p 


12/16 


LC, Ens, IVCL, Tw, Bis 


8 




RLA+VB 


5/16 


LC, Ens 


20 




RLA+NB 


5/16 


LC, Ens 


20 




VA+RLp 


16/16 


LC, IVCL, CR, LD, St, BLs, 
SL 


5 




NA+RLJ3 


16/16 


LC, IVCL, CR, St, BLs, SL 


7 


Tomato 


RLA 


0/16 








RLA + RL(3 


0/12 


_ 


_ 




RLA+VB 


0/12 


_ 


_ 




RLA+NB 


0/12 


_ 


_ 




VA+RLP 


12/12 


SLC, CR, LD, St, BLs, SL, 
LPk 


8 




NA+RLP 


7/12 


LC, CR, YL, Tw 


18 


N. 

benthamiana 


CRA 


16/16 


LC, Ens, Tw 


7 




CRA+CRP 


16/16 


LC, Ens, IVCL, Bis, St, SL 


6 




CRA+VB 


16/16 


LC, Ens, Tw 


7 




CRA+NB 


16/16 


LC, Ens, Tw 


7 




CRA+ RLp 


16/16 


LC, Ens, IVCL, Bis, St, SL 


6 




VA+LKp 


I 0/ I 0 


I C \\IC\ Q\ c T\a; 

LL, IVLL, bLS, I W 


o 




NA+CRP 


16/16 


LC, MY, VC, Bis, 


6 




CRA+VA+CRP 


16/16 


LC, IVCL, CR, Ens, St, BLs, 

T\a; 

I w 


5 




CRA+NA+CRP 


16/16 


LC, Bis, IVCL, Ens, Tw 


5 




CRA+VA+VB 


16/16 


LC, IVCL, CR, Ens, Chls, 

PI r CI T\a/ 

bLS, oL, I W 


5 




CD A i M A i MD 

LKA+NA+Nb 


I 0/ I 0 


1 r VI \/D C+- KA+ Dlr 

LL, YL, Vb, ot, IVIt, blS, 

IVCL, Ens 


5 




CRA+VA+VB+p 


16/16 


LC, IVCL, Ens, LD, St, Bis, 
SL, Tw 


5 




CRA+NA+NB+p 


16/16 


SL, YL, St, LC, Mt, Tw 


5 




VA+VB 


10/10 


LC, IVCL, CR, LD, Chls, Bis, 
SL 


6 




MA i MD 

NA+Nd 


10/10 


C 1 \/l 1 f~ N 

SL, YL, LC, Mt 


5 


Tomato 


CRA 


10/12 


LC, LR 


21 




CRA+CRp 


1 2/1 2 


LC, LR, VC, SL 


18 




CRA+VB 


10/12 


LC, LR 


21 




r~D A i M D 

CKA+Nd 


9/1 2 


LC, LR 


21 




CRA+ RLP 


12/12 


LC, LR, VC, SL 


19 




NA+CRP 


0/12 








VA+CRP 


10/12 


LC, MY 


9 




CRA+VA+CRP 


12/12 


LC, LR, VC, SL, MY 


9 




CRA+NA+CRP 


12/12 


LC, LR, VC, SL 


16 




CRA+VA+VB 


12/12 


LC, LR, CR, YL 


7 




CRA+NA+NB 


12/12 


LC, LR, VC, SL, YL 


8 
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Table 4 Constructs, their infectivity and symptoms induced on agroinoculated plants (Continued) 

CRA+VA+VB+P 1 2/1 2 LC, LR, CR, YL, St 7 

CRA+NA+NB+P 1 2/1 2 LC, LR, VC, SL, YL 8 

VA+VB 10/10 LC, CR, LR, YL 7 

NA+NB 10/10 LC, CR, YL 8 

LC leaf curling; BV big veins; MY mild yellowing; YL yellow leaf lamina; St stunting; LD leaf distortion; VB vein banding; MO severe mosaic; LR leaf rolling; Bis 
blistering; SL small leaves; CR leaf crinkling; IVCL interveinal chlorosis; Chls chlorotic spots; MiMO mild mosaic; Mt mottling; SYL severe yellow leaf lamina; Tw 
twisting of petiole VC vein clearing, LPk Leaf puckering, Ens enation 



[2]. Betasatellites associated with either RaLCV-[IN: 
Var:05] or CYVMV-[IN:Pat:Rad:08] was > 78% identical 
to TbLCB and CroYVMB, respectively, indicating they 
are isolates of previously described begomovirus asso- 
ciated betasatellite species [30]. Hence, betasatellites are 
designated as TbLCB-IN[IN:Var:Rad] and CroYVMB- 
[IN:Pat:Rad:08], respectively, which is also in concurrent 
with results of the phylogenetic analysis (Figure 2). 

A current consensus regarding the extent of begomo- 
virus diversity holds that high frequency of recombina- 
tion, pseudorecombination and mutation is leading to 
emergence of highly pathogenic viruses which in turn 
cause economically important diseases [31-33]. The 
begomovirus causing RaLCD in the Pakistan can infect 
cotton and is strain of CLCuV [23]. However, in the 
present investigation, radish-infecting begomoviruses 
could not infect and produce symptom on cotton till 60 
dpi. This is the first demonstration of Koch's postulates 
of radish-infecting begomoviruses, one of which is pre- 
dominantly associated with infection on weed host, Cro- 
ton. In the context of begomovirus emergence, it has 
been suggested that indigenous viruses infecting weed 
and wild hosts had been transferred to the new host, 
which facilitates emergence of novel species through 
recombination and/or pseudorecombination events 
[31-33]. Several new begomovirus species isolated and 
identified have arisen by recombination [28,34]. 

Recombination can provide selective advantage in the 
evolution of viruses at species, genera and family levels 
and has been documented in both animal and plant 
viruses with RNA or DNA genomes [35-37]. In this 
study, recombination was detected among full-length 
DNA-A component of one of these newly identified rad- 
ish-infecting begomovirus species and it's near relatives 
by RDP analysis. The pattern of recombination and posi- 
tion of breakpoints is very convincing and supports ear- 
lier study about pattern of recombination in 
geminiviruses [32,38]. Breakpoint was detected in the vir- 
ion sense ORF AV1/AV2, in RaLCV-[IN:Var:05] genome 
and for this region it shared maximum identity with PaL- 
CuV. ORFs AC2 and AC3 of RaLCV-[IN:Var:05] also 
shares maximum identity for this region with PaLCuV. In 
the rep-coding region, recombination was found with a 
weed infecting begomovirus most similar to EuLCV. This 



region is commonly exchanged among begomoviruses 
[38] and may provide a selective advantage, such as 
enhanced replication or suppression of host defense. It is 
apparent that RaLCV-[IN:Var:05] is a chimeric molecule 
and has arisen by inter-specific recombination among 
EuLCV, PaLCuV and ToLCBDV. This result was also 
confirmed through Simplot analysis. Further, placement 
of RaLCV-[IN:Var:05] in a distinct phylogenetic tree with 
isolates of EuLCV, PaLCuV and ToLCBDV, monopartite 
begomoviruses from Asia supported this observation. 
These findings together indicate that recombination 
among begomoviruses is an important driving force in 
evolution of plant infecting begomoviruses with extend- 
ing host range. An increasing number of monopartite 
begomoviruses has been reported in India in recent 
years, most of them are associated with DNA betasatel- 
lites molecules with evidences that some of them may 
have not evolved in the recent past [34]. 

Agroinoculation studies demonstrated that these 
monopartite viruses (RLA and CRA) can replicate 
autonomously and move systemically in their natural 
host (radish), and induce mild leaf curl symptoms. How- 
ever, in the presence of DNA betasatellite, RLA or CRA 
induces these symptoms, which appear similar to those 
observed in the fields in Varanasi and Pataudi, India. 
Therefore, the RaLCD disease complex is caused by cog- 
nate association of DNA-A and DNA betasatellite. Pre- 
requisite association of a begomovirus and satellite 
DNA-P has also been found to be necessary for typical 
disease development [14,15,34]. However, association of 
begomovirus and betasatelite complex with RaLCD has 
been demonstrated for the first time. Inability of RLA or 
CRA to transreplicate DNA-B molecule from either 
ToLCGV or ToLCNDV is attributed due to difference 
in iteron sequences and Rep binding sites observed. 
However, a number of recent studies have demonstrated 
that Rep-mediated binding can be more relaxed than 
thought previously [29,39]. Increase in symptom pattern 
and severity associated with betasatellites could be due 
to a direct effect of the PCI protein on plant develop- 
ment [40,41]. Trans-complementation of RLp and CRp 
with CRA and RLA, respectively, revealed the promiscu- 
ous nature of betasatellite replication. This was further 
evidenced from the finding that two tomato-infecting 



Singh et al. Virology Journal 2012, 9:43 
http://www.virologyj.eom/content/9/1/43 



Page 13 of 18 



Nicotiana benthamiana 



Tomato 



ca. „ + +++ + + 

0 3?z>z>z>zx 



Vfif 



8" 5 S f 



8 8 



5? ^ 



a) 



s 




J88 


8 8 i iisi 








'•Mf. 





100 88 76 80 78 74 86 88 68 141 103 



100 100 100 97 86 95 100 99 93 140 100 



b) 



86 89 49 58 53 42 



c) 



31 102 115 49 100 



i i i i .u^iSLri ■ ■ ■ is 

100 101 100 





97 98 80 78 70 



- 6 



Probe 

■CYVMV-A 

HCYVMV-p 
HoLCNDV-A 



68 98 



100 



-♦ToLCNDV-B 



101 100 100 




HToLCGV-A 



40 20 100 100 100 



40 20 110 90 100 - 




<*ToLCGV-B 



Figure 5 Relative level of DNA accumulation following reassortment and mixed infection of begomoviruses. Southern blot analysis 

showing relative accumulation of viral DNA and at the bottom of gel the numerical value is showing in percentage of viral DNA accumulation 

in inoculated tomato (right panel) and N. benthamiana (left panel). 
^ ) 



begomovirus species (ToLCGV[IN:Var:01] and 
ToLCNDV[IN:Svr:93]) also served as helper viruses for 
RLp and CRp in N. benthamiana [42-45]. Our results 
also show that betasatellites complement the functions 
of the DNA-B component of tomato-infecting 



begomoviruses (ToLCGV and ToLCNDV), by mobiliz- 
ing DNA-A from sites of inoculation to the distal tis- 
sues. These results fortify the idea that the PCI protein 
can functionally substitute the role of DNA-B move- 
ment proteins [46]. 
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Mixed infection of begomoviruses is common in 
Indian sub-continent and situation is further worsened 
owing to ability of whiteflies to transmit more than one 
begomovirus species under field conditions [47,48]. 
Since, Croton yellow vein mosaic virus and betasatellite 
DNA, isolated from RaLCD infected plant could infect 
tomato, we further extended study to test interaction of 
CRA and CRp with predominant tomato-infecting bego- 
moviruses species (ToLCNDV and ToLCGV) from 
India. Plants (N. benthamiana and tomato) inoculated 
with mixtures of either combinations of CRA + VA + 
CRp, CRA + NA + CRP, CRA + VA + VB, CRA + NA 
+ NB, CRA + VA + VB + CRp and CRA + NA + NB + 
CRp, produced severe symptoms compared with plants 
infected with isolates of only one begomovirus species. 
Viral DNA accumulation (CRA) was also slightly 
enhanced in plants inoculated with CRA + VA + VB + 
CRP or CRA + NA + NB + CRp inoculated plants com- 
pared to plants inoculated with CRA + CRp. Interest- 
ingly, beta satellite accumulation was decreased (upto 
60%) in the presence of DNA-B component (either NB 
or VB) in the above combinations tested. This is in con- 
current with earlier observation [49]. No significant 
increase of DNA accumulation of NA or VA was 
observed while increase of CRA in dual infected plants 
indicates asymmetric synergism (Figure 5). Asymmetric 
synergism has earlier been reported among two tomato- 
infecting begomoviruses [29]. 

Conclusion 

RaLCD in India is caused by association of two monopar- 
tite begomoviruses and promiscuous betasatellites. This is 
also the first experimental demonstration of Koch's postu- 
late for begomoviruses associated with leaf curl disease of 
radish. The aetiology of this disease and the capacity of 
the DNA- betasatellite to be replicated by other begomo- 
virus species increase the possibility of recombination and 
reassortment events. This could lead to evolution of new 
recombinant viruses or begomovirus complexes with dif- 
ferent biological properties. It is also possible that 
exchange of betasatellites (RaLCV associated with TbLCB 
satellite DNA-P) could extend the virus host range thereby 
emergence of new diseases in cultivated crop plants (for 
example tomato, radish). The persistence of promiscuous 
betasatellites, such as RLp and CRp in radish or other 
hosts like weeds (e.g. Croton spp.) could facilitate the 
emergence of highly pathogenic begomovirus/betasatellite 
complexes that may induce more severe diseases by over- 
coming natural resistance of hitherto non-host. 

Methods 

Sample collection, grafting and DNA extraction 

Surveys of radish-growing areas of Varanasi, Uttar Pra- 
desh and Pataudi, Haryana in India were conducted and 



infected radish plants showing typical virus infection 
symptoms (stunted and distorted growth and leaf curl) 
were collected. The leaves alongwith petioles of at least 
1-1/2 in. long were cut from the symptomatic plants 
and were used as scion. Infected scions were grafted 
onto the healthy radish leaf petiole, known as stock. 
Grafted plants were kept within an insect-proof glass 
house. Total nucleic acids were extracted from collected 
samples [50]. 

Viral DNA detection and cloning 

General detection of begomovirus association was car- 
ried out by Southern hybridization using coat protein 
probe [29] and by PCR with degenerate begomovirus 
specific primers [25]. PCR amplified DNA fragments 
were cloned into the pTZR57R/T vector (Fermentas Inc, 
USA) and sequenced. To facilitate full-length cloning of 
RaLCV- Varanasi genomic component (DNA- A), unique 
restriction site Nde I was determined, around which pri- 
mer pairs RINdelFP (S'CATATGTGAGCCGTGTTGS') 
and RINdeIRP (S'CATATGGGCTGTCGAAGTS') were 
designed. 

For cloning of RaLCV- Pataudi, Haryana samples, viral 
DNA were amplified using Phi 29 DNA polymerase 
(Temphilphi; GE Healthcare) according to manufac- 
turer's instruction. Briefly, the amplified concatemers 
were digested with restriction endonucleases to yield 
presumed monomeric virus components of -2.8 kb. 
Potentially, full-length linearized viral genomes (at Bam 
HI site) were cloned into pKS + vector. 

DNA-B degenerate primers [25] were used to amplify 
DNA-B component. The inability to amplify DNA-B 
using specific primer led us to design new set of primers 
corresponding to the common region shared by DNA-A 
and DNA-B followed by restriction endonuclease diges- 
tion of purified super-coiled DNA to identify two popu- 
lations of molecules. Presence of betasatellite molecules 
was detected by PCR with primers beta01/02 [26]. 
Amplified products (-1.3 kb) were cloned into 
pTZR57R/T vector (Fermentas Inc, USA). 

Sequence comparisons 

Putative full-length clones were purified using a Qiagen 
plasmid miniprep kit (Qiagen) and their sequences were 
determined commercially (Chromous Biotech, India). 
The sequence results were analyzed by comparing the 
sequence to other begomovirus sequences in GenBank 
using BLAST, followed by analysis using MEGA pro- 
gram version 4 (http://www.megasoftware.net/mega4/ 
features.html) [51] with default parameters. The full- 
length genome sequences of representative begomo- 
viruses and betasatellites were aligned using ClustalW 
and dendogram was generated using 1000 bootstrap 
repetitions. Standard abbreviations and GenBank 
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accession numbers for DNA-A [2] and DNA betasatel- 
lites [30] were used. 

Recombination analysis 

Detection of potential recombinant sequences, identifi- 
cation of likely parental sequences, and localization of 
recombination breakpoints were carried out with the 
recombination detection program (RDP version 3) [52]. 
Default RDP settings were used throughout with a P- 
value cutoff of 0.01 and the standard Bonferroni correc- 
tion. For bootscan analysis, 200 replicates with a 95% 
cut-off were taken, and for GENECONV analysis, the g- 
scale parameter was set to 1. In addition, both the viral 
genome and DNA betasatellites sequences were analysed 
using Simplot program, version 3.2 (http://sray.med.som. 
jhmi.edu/SCRoftware/Simplot), with a sliding window of 
200 nucleotides moving in 20-nucleotide steps [53-55]. 

Generation of constructs for infectivity studies 

To check infectivity, begomovirus and betasatellite were 
cloned as partial tandem repeat in pCAMBIA2301. For 
RaLCV-Var, -1.6 kb Bam HI (152)-£co RI (1789) frag- 
ment containing the intergenic region (IR) was cloned 
to generate a 0.6 mer (pRLAV0.6). The full-length 
monomer was cloned into the Bam HI digested 
pRLAV0.6 to generate a 1.6 mer tandem repeat called 
pRLAV1.6 (referred to here as RLA). For RaLCV-Patudi, 
a 2.4 kb Bam HI-Hind III fragment containing the IR 
was cloned to generate -0.8 mer (pCRAP0.8). The full- 
length monomer was cloned into the Bam HI digested 
pCRAP0.8 to generate 1.8 mer tandem repeat called 
pCRAP1.8 (referred to here as CRA). 

For the RaLCB-Var, Pst I (56)-Kpn I (1298) fragment 
was cloned into pCAMBIA2301, followed by ligation of 
full-length Kpn I fragment to produce partial tandem 
repeat called pRLpVl (referred to here as RLp). For the 
RaLCB-Pat, Kpn I (1292)-£co RI (330) 395 bp fragment 
was released from RaLCB-Pat DNA-P and was cloned 
into pCAMBIA2301, followed by ligation of full-length 
Kpn I fragment to produce 1.3 mer tandem repeat called 
pCRpP1.3 (referred to here as CRP). Insert orientation 
within partial tandem repeats were confirmed by restric- 
tion digestion with appropriate enzymes. 

Infectivity and trans-complementation of cloned DNAs 

Recombinant plasmids (RLA, RLp, CRA and CRp) and 
the 'empty' vector (used as negative control) were intro- 
duced into Agrobacterium tumefaciens strain EHA105 
by transformation. Agro-inoculation was performed in 
radish (cv Japanese White), Nicotiana benthamiana and 
tomato plants. Young seedlings (3-4 leaf stage) were 
inoculated viral constructs either alone or in combina- 
tion by needle puncture method. 



Trans-complementation with genomic components of 
tomato-infecting begomoviruses {Tomato leaf curl New 
Delhi virus [ToLCNDV] and Tomato leaf curl Gujarat 
virus [ToLCGV]) [32] was also studied. Agrobacterium 
cultures containing respective constructs (Table 4) were 
mixed in equal concentration and inoculated into test 
plants. The inoculated plants were placed in a controlled 
growth chamber with 16 h day light at 28°C/70% relative 
humidity. All the experiments were repeated at least for 
four times and symptom progression was recorded till 
45 dpi. 

Detection of viral DNA 

Total genomic DNA was extracted from leaf tissue at 30 
dpi following Dellaporta method [50]. Viral and satellite 
DNAs was detected by PCR with the specific primers 
for RLA (RLF1.3 5XTGGGCTTACCCATAGAGTGG3' 
& RLR1.3 S'CAGGGAAGACAATGTGGGCCTS'), CRA 
(CRF0.9 S'ATGGGTCTCTGCATATCCATGCCCTCS' 
& CRR0.9 5TCAACTCGTCGACGCCTGATCCC 
TTTC3'), (RLPF 5'ATGACGATCAAATACAAAAAC- 
CAGAAAG3' RLpR 5TTATACAGATGA ACGCGTA- 
TACACATCG3') CRp (CrpF 5'ATGACGATCATA 
TATCAGAATGAGAC & CrpR 5'TTACACATTTACA- 
TATTTAGACACATC3'). In addition, specific primers 
were used for detection of genomic components of 
ToLCGV and ToLCNDV following Chakraborty et al. 
[29]. 

For southern hybridization analysis, total genomic 
DNA (8 (ig) was separated on agarose gel (1%) and trans- 
ferred to Hybond-N+ membrane (Amersham GE Health- 
care, UK), following standard procedures [56]. Viral 
DNA was detected by hybridizing blots separately using 
radiolabeled probes of DNA-A and DNA- betasatellites 
specific to either species. For RLA, a BseRI fragment (nt 
2097-2616) and for RLp a EcoRI fragment (nt 299-506) 
was used as the probe. For CRA, a Pst I-Xba I (nt 1455- 
1884) fragment and for CRp a EcoRI-Nde I fragment (nt 
330-1041) were used. Specific probes for detection of 
DNA-A and DNA-B of ToLCNDV and ToLCGV were 
prepared according to Chakraborty et al [29] . DNA frag- 
ments were labelled with [a- 32 P]dCTP by random oligo- 
nucleotide-primed synthesis [57]. Viral DNA levels were 
quantified using a Phosphorlmager (Fuji Film). 
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